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I : Abstract. A model interpreting variations of EUV brightness upward propagating in 

solar coronal loops as slow magnetoacoustic waves is developed. A loop is considered to 

^~^ ■ have a non-zero plane inclination angle and offset of circular loop centre from the base- 

On line. The model also incorporates effects of dissipation and gravitational stratification. A 

I^ ■ . 

l/~j , linear evolutionary equation is derived and applied to investigation of slow wave dynam- 

r^ . 

r-^ ics. Both the non-zero plane inclination and the offset decrease the growth of the wave 

^.^ ' amplitude due to stratification. It is shown that wide spectrum slow magnetoacoustic 

(-H [ waves, consistent with currently available observations in the low frequency part of the 

"^i: spectrum, can provide the rate of heat deposition sufficient to heat the loop. In this 

5^ scenario, the heat would be deposited near the loop footpoints which agrees with the 
-*— > 
^ current observational data. 

> 

• I— I 
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1. Introduction 

Identification of possible mechanisms of heating of solar coronal plasma and acceler- 
ation of solar wind is among the unsolved problems of solar physics. In the last two 
decades there has been intensive theoretical work undertaken to advocate the idea of 
coronal plasma heating by magnetohydrodynamic (MHD) waves. In particular, it has 
been suggested that coronal loops could be heated by dissipation of high frequency slow 
magnetoacoustic waves (Porter, Klimchuk & Sturrock 1994a, b; Laing & Edwin 1995). 
However, up until recent observational advances, these studies have been poorly moti- 
vated because of the absence of observational evidence of the MHD wave activity in the 
solar corona. 
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Recent observations performed with SOHO/EIT and TRACE imaging telescopes have 
established: 

1. Existence of upward propagating perturbations of EUV brightness in long coronal 
loops at the bandpasses EIT 195A (Berghmans & Clette 1999) and TRACE 17lA (Bergh- 
mans et al. 1999, Dc Moortcl, Ireland & Walsh 2000), and EIT 195A and TRACE 171 A 
simultaneously (Robbrecht et al. 2001). Probably, the same phenomenon was observed 
by Nightingale, Aschwanden & Hurlburt (1999) in TRACE 17lA . 

2. The projected velocities of these plasma density perturbations are about 65-165 
km/s, which is below the sound speed corresponding to the temperature of the loop. 

3. The amplitude is usually 1-3% of the background density, reaching 5-6% in some 
observations. 

4. The propagating disturbances are quasi-periodic with the periods of 3-20 min. 
There is no a well-defined correlation between the periods and physical parameters of 
the loops detected. 

5. The disturbances are detected in the loop legs near footpoints and always propagate 
upwards. There are no observations of waves propagating downwards. 

Two rival interpretations of the propagating disturbances have been suggested. In the 
first (Reale, Peres & Serio 2000), the disturbances are considered as field-aligned flows 
of matter, generated by the siphon effect. In the second (Nakariakov et al. 2000), the 
propagating disturbances are associated with slow magnetoacoustic waves. The wave- 
based theory is similar to the model developed for the interpretation of compressive 
perturbations observed in polar plumes (Ofman, Nakariakov & Deforest 1999, Ofman, 
Nakariakov & Sehgal 2000) in terms of slow magnetoacoustic waves. Both theories are 
consistent with the observations and with the MHD theory, but the flow-based model has 
difficulties with explanation of the periodicity and the symmetry of the front and back 
edges of the disturbance. Also, the observations do not show a change (acceleration) in 
the propagating disturbance speed, predicted by the flow-based theory. In contrast, the 
wave-based theory explains all the observational flndings and, in particular, justifies the 
absence of downward propagating disturbances. 

In the wave-based theory (Nakariakov at al. 2000), it has been established that the 
main physical factors, which are important in the slow magnetoacoustic wave dynamics in 
long coronal loops, include: (i) gravitational stratification, which is responsible for wave 
amplitude amplification; (ii) dissipation, which leads to the conversion of wave kinetic 
energy into heat in the large wave numbers part of the spectrum, and (iii) nonlinearity, 
which generally yields wave front steepening by generation of higher harmonics, leading 
to enhanced dissipation. An evolutionary equation for the slow magnetoacoustic wave 
amplitude has been derived, which includes these three effects. However, in order to catch 
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the physical essence of the problem without further complications, simple geometry of 
the loop, namely, zero inclination angle and zero offset of the centre from the solar surface 
has been assumed. 

In contrast, observations show that the plane of a loop is not necessary perpendic- 
ular to the Sun's surface (a non-zero inclination) and the centre of the loop circle is 
usually situated below the surface (a non-zero offset) (e.g. Aschwanden et al. 1999). Con- 
sequently, a more elaborate theory of the slow magnetoacoustic waves in coronal loops, 
which includes the effects of the non-zero inclination and the non-zero offset, is required. 

Soft X-ray and EUV coronal telescopes provide the possibility to determine the distri- 
bution of density, temperature and heat deposition along a loop. Priest et al. (2000), using 
Yohkoh/SXT observations concluded that hot, soft X-ray loops were heated uniformly. 
Aschwanden et al. (1999), Aschwanden, Nightingale & Alexander (2000) and Aschwan- 
den, Schrijver & Alexander (2001) analyzing SOHO/EIT and TRACE EUV observations 
argued that cooler loops were heated non-uniformly, near the loop footpoints. The spe- 
cific mechanism for the energy deposition remains unknown and could, in principle, be 
associated with dissipation of slow magnetoacoustic waves. 

De Moortel et al. (2000) have estimated the energy of slow waves observed to be 
insufficient to heat a loop. This result was based upon the assumption of the wave to be 
a single harmonic with the period of several minutes. But, currently available resolution 
of the telescopes does not allow us to detect periodic motions with the periods shorter 
than a few minutes, even if they are presented in the spectrum. On the other hand, the 
absence of well-defined period in the loop waves suggests that their spectrum is probably 
continuous and not a single harmonic. We would like to point out that the main difficulty 
with the theory of coronal loop heating by dissipation of slow magnetoacoustic waves has 
been connected with the absence of the observational evidence of these waves in the 
corona. It has been believed that the waves somehow generated below the temperature 
minimum, are not able to penetrate through this layer. Now, the observational evidence 
of the long period longitudinal waves has been provided. The long period waves can have 
unresolved shorter period components in their spectrum. According to the theory, these 
short period spectral components dissipate very quickly, and deposit the energy in the 
near footpoint part of the loop legs. 

In this paper we develop further the model of Nakariakov et al. (2000) for one- 
dimensional slow magnetoacoustic waves in stratified coronal loops, taking into account 
the effects of the non-zero loop plane inclination angle and the semi-circular loop off- 
set, and derive a generalized evolutionary equation for slow magnetoacoustic waves in 
coronal loops. It is shown that a coronal loop acts as a dissipative filter on the slow mag- 
netoacoustic waves and that the short wavelength part of the wave spectrum dissipates 
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Fig. 1. The sketch of the model considered. A coronal loop is considered as a straight 
magnetic field line with density and gravitational acceleration varying along the axis of 
the cylinder. 



in the nearest vicinity of the loop footpoints. By calculating the heat deposition rate, 
provided by the dissipation of the slow waves and comparing it with the observationally 
determined heat deposition rate in EUV loops, we show that sufficiently wide spectrum 
slow waves can deposit the required amount of energy in the lower parts of the loops, 
which is consistent with observations. 

This paper is organized as follows: we formulated our generalized model in Section 2; 
Section 3 deals with the comparison of our model's predictions with the observations; 
and, finally, in Section 4 we close with a discussion of our main results. 



2. The model 



We consider a semi-circular loop of the curvature radius R^. The loop cross-section is 
taken to be constant. Validity of the latter assumption has been corroborated by Klim- 
chuk et al. (1992) and by Aschwanden et al. (1999) for soft X-ray and EUV loops, 
respectively. The loop is filled with a gravitationally stratified magnetized plasma of con- 
stant temperature. Here, we generalize Nakariakov et al.'s (2000) work by incorporating 
effects of loop plane inclination and offset of circular loop centre from the baseline. This 
is necessary since actually observed loops have the non-zero (measured from the normal 
to the solar surface) inclination angle a and non-zero offset of circular loop centre from 
baseline, Zq, i.e. distance from the circle's centre to the solar surface. The sketch of the 
model geometry is shown in Figure 1. 
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Fig. 2. Evolution of amplitude of a harmonic slow magnetoacoustic wave propagating 
along a coronal loop of the radius 93 Mm and zero off-set for different loop plane incli- 
nation angles. The wave period is 300 s. The medium is dissipationless, fj — 10^^. The 
solid curve corresponds to the inclination angle a = 0, and long and short dashed lines 
represent a = 7r/4 and a = 7r/3, respectively. 



The gravitational acceleration projected on the loop tangential at any point which is 
located at a distance s measured from the footpoint (s — 0) along the loop, is 
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where G is the gravitational constant, i?^ is the loop radius, Rq and Mq are solar radius 
and mass, respectively. 
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Fig. 3. Evolution of amplitude of a harmonic slow magnetoacoustic wave propagating 
along a coronal loop of the radius 93 Mm and zero inclination angle. The wave period 
is 300 s. The medium is dissipationless, 77 — 10"^. The solid curve corresponds to the 
loop with the off-set Zo = 0, and long and short dashed lines represent Zo = Rl/^ and 
Zo = Rl/^ respectively. 



The use of the hydrostatic equilibrium equations and the isothermal equation of state, 
allows us to write the stationary density profile along the loop as 



Po(s) = po(0) exp ( -— / g{s')ds' 



(2) 



where 7 = 5/3 is the adiabatic index and Cs is the speed of sound. The temperature 
T was assumed to be 1 MK, giving the sound speed Cs of 152 km/s. An effective scale 
height can be determined as H{s) = C'^/jg{s). 

Inserting generahzed expression (|^) for g{s) into (||) we obtain 



Po(s) =po(0)exp 



75(0) 
C? 
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Fig. 4. Evolution of the slow wave amplitude in the average EUV loop (Aschwanden et 

o 

al. 1999) with the parameters: Rl ~ 93 Mm, Zq = 62 Mm, a = 7 . The normalized 
coefficient of dissipation fj is 4.0 x 10^'*. The solid curve corresponds to the wave period 
of 600 s, while long dashed, short dashed and dotted curves correspond to 300 s, 120 s 
and 60 s respectively. 



RL^l-{Zo/RL)^sm{s/RL) - Zo{l - cos(s/i?i)) 



cos a 



1 + cos a 



(3) 



{Rl/Rq)^/! - (Zo/Rl)^ sin(s/i?L) - {Zo/Rq){1 - cos{s/Rl)) 

Note, that when a — > 0, expressions (P and (^) reduce to the case investigated by 
Nakariakov et al. (2000). 

For the following, we introduce dimensionless variables, Rl = Rl/ Rq, Zq = Zq/Rq, 
p = p/ pvi, H = H/Rq, V = V/Cs and s — s/Rq. Assuming that the wavelength is 
much less than both the scale height and the dissipation length, and applying the multi- 
scale expansion method described in detail in Nakariakov et al. (2000), we obtain the 
evolutionary equation for slow magnetoacoustic waves, 
dV 1 ^-, fj d^V 



ds 2H{s)^ 2p{s) de 



= 0, 



(4) 
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where ^ = s — Cgt/ Rq is the running coordinate and 



1 
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3 7^7 



(5) 

is the normalized coefficient of dissipation connected with the coefficient of thermal con- 
duction Km and the volume viscosity rjQ (see the discussion in Nakariakov et al. 2000), 
the rest of the notations are standard. 

Eq. fl) together with ffl) and (0) generalizes the equation derived by Nakariakov et 
al. (2000) to the case of non-zero loop plane inclination and the offset, Zq. But, it does 
not contain the nonlinear term. This term is neglected because it was shown that, for 
the amplitudes observed, the effect of the nonlinearity is insignificant. Neglecting the 
nonlinearity allows us to investigate solutions of the evolutionary equation in all the 
details analytically. 

Considering a harmonic wave, oc exp{ik£,), we solve Eq. dJ) for the amplitude of the 
fc-th spectral component, 

V{s, k) = ViO,k) exp [01 (s) - fc202 (s)] , (6) 

where 

"^'^ ' Jo 2H{s') 



( 75(0) 



K'^C^sRq 



Rl^I--{Zo/RlYsHs/Rl) - Zo{l - cos{-s/Rl)) 



cos a 



/ 1 -f cos a 



Rl^I-{Zo/RlY^Ms/Rl) - Zo{l - cos{s/Rl)) 



(7) 



M-^)= / TT^Jrs', (8) 

Jo Ms) 

and k is the dimensionless wave number, k — 2ttRq/X — 2'iiRq/ {CsP) with A and P 
being the wavelength and wave period respectively. Here, we used the fact that short 
wavelength slow magnetoacoustic waves are practically dispersionless and described by 
the dispersion relation lo « Cgk. Consequently, their spatial and temporal spectra are 
similar. 

The parametric studies of solution (o) are shown in Figures 2 and 3. Fig. 2 demon- 
strates the behaviour of a harmonic wave in identical loops with different values of the 
inclination angle a. Naturally, the increase of the inclination angle decreases the value of 
the maximum amplitude the wave can reach. Indeed, larger inclination angle means less 
stratification of the plasma density in the loop, therefore lower maximal values. Fig. 3 
shows the evolution of the wave amplitude along the loop for different offsets of a circu- 
lar loop centre from the baseline, Z^- The increase in the loop centre offset leads to the 
decrease of the loop length and also decreases the angle between the loop leg and the 
surface of the Sun at the footpoint. The length of the loop is L = 2Rl arccos(Zo/i?L)- 
Obviously, the increase in Zg reduces the effective stratification of the loop density and 
this results in lower maximal amplitudes reached by the waves. 
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3. Evolution of the spectrum 

Eq. (|y) describes evolution of the spectral component k. In the exponent, the first term 
is positive in the ascending stage of the wave propagation, becomes zero when the wave 
reaches the loop apex and then becomes negative. The second term is always negative. 
Consequently, in the ascending stage, the wave evolution is described by the competition 
of these two terms describing effects of stratification and dissipation, respectively. Im- 
portantly, these two mechanisms operate in different spectral domains. Mathematically, 
Eq. (jm is similar to an evolutionary equation of a dissipative filter in the theory of radio 
circuits. Therefore, a coronal loop acts on slow magnetoacoustic waves as a dissipative 
filter. If we knew the initial spectrum of the wave, the equation would allow us to deter- 
mine the spectrum at each given height s of the loop. Unfortunately, the initial spectrum, 
especially its short wavelength part, cannot be determined with the current resolution of 
EUV telescopes. 

Evolution of the spectrum is accompanied by the dissipation of spectral components. 
Different spectral components have different effective scale heights of the dissipation. 
Assuming the initial wave spectrum and following the spectrum evolution along the 
loop, we can determine the energy deposition rate at each given height s. Indeed, this 
function can be tested against the observations. 

In order to compare predictions of our model with the observational data, we have 
investigated the behaviour of different slow magnetoacoustic wave spectral components in 
the average EUV loop determined by Aschwanden ct al. (1999), with Rl = 93 Mm, Zq = 

o 

62 Mm, a = 7 , corresponding to the loop length of about L — 156 Mm. The evolution 
of different spectral components is shown in Figure 4. The dissipative coefficient fj has 
the same value as was estimated in (Nakariakov et al. 2000). The spectral components 
corresponding to shorter wavelengths (i.e. shorter period) dissipate faster and closer to 
the loop footpoint. 



The initial spectrum of the slow waves, V{0, k) = \J f{k) is unknown, but it should 
satisfy certain criteria. According to observations, the spectrum should have a maximum 
in the vicinity of 180-420 s. The velocity amplitude of the maximum is about 1-3% 
of the speed of sound. Perhaps, there are also maxima in the shorter period part of 
the spectrum, but they have not been resolved observationally. So, we assume that the 
spectrum gradually declines with the wave period. In order to avoid singularities in the 
integration of the dissipated wave energy over the whole spectrum k — [0, -foo), we 
assume that /(fc) -^ when k — > -|-oo and fc — > 0. Thus, we choose the spectrum /(fc) 
with a single maximum at A:* . A function which satisfies these requirements is 

/(fc)=n"/2(^)^ '' ki^a + kA (n>3). (9) 
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Fig. 5. Possible initial spectra f{k) of slow magnetoacoustic waves. The solid curve cor- 
responds to n = 3, and long and short dashed curves to n = 5 and 7, respectively. The 
diniensionless wave number k is connected with the wavelength X as k ~ 2ttRq/X. 

The spectra with n ~ 3, 5, 7 are shown in Figure 5. For any n, the spectrum /(fc) attains 
its maximal value at A:* = \/a/2. For a — 1.84 x 10^, fc, « 96, which corresponds to 
a wave period of 300 s, and n = 3, the actual spectral peak {f{k) > 0.69) is situated 
between the wave numbers corresponding to 1.5 and 12 min periods, which is consistent 
with observations. Note that f{k) ex fc^*^"^^^ asymptotically as fc ^ +oo. The energy 
contained in the spectrum is determined by the spectral index n. Higher values of the 
spectral index correspond to narrower spectra. 

The wave energy flux dissipated into heat in the loop segment from the footpoint to 
the distance s is defined as 

r + oD 

Fd{s) - / [FokiO, k) - Fokis, k)] dk, (10) 

Jo 

where 



FDk{s,k) ^ -po{s)V^iO,k)A^C^^cxp[2Ms)-'2k''Ms)] 



(11) 
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and A is the amplitude of the strongest harmonic in the spectrum, measured in the units 
of the sound speed Cg, and Eq. (p) was used. To evaluate the integral in Eq. dlO), we use 
the spectrum (0) and expressions (M) and (H). For the odd integer n, the integrals can be 
evaluated analytically, in particular 

+00 



FoTis 



k)dk = ^poA^C^^{1.5V3a)(a~^/^ - ^/2^T(t)2{s)e^'''^^''>'' 



(12) 



-t-Ctt 



erf iy'2a4>2{s) 
F5f (s,fc)dfc-ipoA2C,3(5^/2(a/2)V3)('(2a)-^/2^20F(2(/.2(s))3/2e2'^^(^)(2-)(i3) 

1 - erf (v/(2a)202(s)) 



2(202(s))(2a)-V2 



and 



r+00 1 / 



1 - erf (V(3a)202(s)) 



2(2(/.2(s))(3a)-3/2 + 4(202(s))2(3a)-i/2 



Here, erf(a;) = (2/y^) J^ exp(— i^)(ii is the error function. In the evaluation of the inte- 
grals, we use that po(s) exp[2(/)i(s)] = 1 and 



I{n ^ 3) 



I{n ^ 5) 



e-^'''/^{x + aYAx = 2/7^- 20i7Ie'^'^(l - erf(^/^)), 



e-^7V(a; + 2a)5da; = - ( (2a)-3/2 + 2V^^3/2gaM(i , 



(15) 
(16) 



erf(^/^)) - 2/i(2a)"i/2 
I{n = 7)= e-^7yi^r^3^da;=— 3(3a)-5/2_4^^5/2gaM(i_ (^7) 

erf(V^)) - 2/i(3a)~^/2 -K V(3a)-i/2 J 

The observationally determined rate of heat deposition at a position s of the loop can 
be approximated as 



Eh{s) = Eho exp{-s/sH), 



(18) 



where Ehq is the heating rate at the footpoint and sh is the heating scale height. The 
specific values of Eho and sh were determined by Aschwanden et al. (1999) for EIT loops 
and by Aschwanden et al. (2000, 2001) for TRACE loops. A physical value identical to 



F£i{s) introduced in this paper, can be constructed by integration of Eq. (18) over s, 

Flj'%s)= f EH{s)ds = EHoSH{l-e-'/'"). (19) 
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Fig. 6. The slow magnetoacoustic wave energy flux dissipated into heat in a loop segment 
from the footpoint to the position s in comparison with the observationally determined 
nonuniform heating of an average loop. The thin solid curve corresponds to the spectral 
index n = 3, the long dashed curve to n = 5 and short dashed curve to n = 7. The 
observational dependence is the thick solid curve with the error bars. 



When s = L/2, at the loop apex, the value of the function Fg'"'(s) is the total energy 
power introduced by Aschwanden et al. (2001, Eq. (12)). 

The comparison of the theoretical dependences of the wave energy flux dissipated 
into heat in a loop segment from the footpoint to the position s given by Eqs. (|l2), (O) 
and (n3) and the observationally determined dependence (n9) is shown in Figure 6. Here, 
we use the parameters of nonuniform heating determined in (Aschwanden et al. 2001), 
Eho — 10^^°^"'^ erg cm^'^ s^-^ and s^ = 12 ± 5 Mm. The normalized dissipation is 
fj — 4.0 X lO"*. The density of plasma near the loop footpoint, which corresponds to the 
concentration at the base of the average loop, no = 1.92 x 10^ cm^'^, is po = 4.1 x 10^^'^ 
g/cm"^. The amplitude A of the wave, corresponding to the peak in the spectrum (0), is 
taken to be ^ = 0.03, i.e., 3% of the sound speed. The observational curve has the error 
bars determined by choosing plus and minus signs in the Eho ^^^d sh, given above. 
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Note that for n = 3 the theoretical curve -Fd(s) fits the median value for the F^''{s) 
very well and that the curves corresponding to much steeper spectra, with n = 5 and 
n — 7 are also almost inside the error range. 



4. Discussion 

Recently observed variations of EUV emission intensity, propagating upward along the 
legs of solar coronal loops, are satisfactory interpreted in terms of the slow magnetoa- 
coustic waves. We generalized the model of the slow wave evolution in a loop, suggested 
by Nakariakov et al. 2000 by taking into account effects of loop plane inclination and 
offset of the semi-circular loop centre from the solar surface. An evolutionary equation 
which governs the wave propagation in presence of gravitational stratification and dissi- 
pation was derived. Both non-zero inclination and offset affect the slow wave evolution, 
reducing the stratification of the loop plasma and, consequently, the amplitude growth 
of the upward propagating waves. 

It was found that coronal loops act as dissipative filters of slow magnctoacoustic 
waves. Different spectral components have different dissipation lengths and the wave 
spectrum evolves with height or with the distance along the loop. Therefore, the amount 
of wave energy deposited in the loop decreases with the height and the wave dissipation 
deposits the energy in the lower parts of the loops. Making certain assumptions about the 
spectrum of the waves near the loop footpoints, we constructed a function describing the 
wave energy fiux dissipated into heat in a loop segment. This function was compared with 
an observationally determined function approximating the nonuniform heat deposition 
requirements. We conclude, that the observed heat deposition into coronal loops can be 
interpreted as dissipation of slow magnctoacoustic waves, provided the wave spectrum is 
sufficiently wide. This conclusion is consistent with previously obtained results by Porter 
et al. (1994a, b), but is based upon recently obtained observational information on the 
actual existence and parameters of slow magnctoacoustic waves in coronal loops and on 
the nonuniform distribution of heat deposition in EUV loops. In the scenario proposed, 
the short wave length components of the slow wave spectrum, with short dissipation 
lengths, are filtered out in the loop legs, depositing the energy and heating the loops near 
the footpoints, while the long wave length spectral components, with longer dissipation 
lengths, propagate further and are detected by EUV telescopes. 

Strictly speaking, realistic spectrum of slow magnctoacoustic waves is restricted by 
the acoustic cutoff wave number from below and by the inverse mean free path length 
from above. However, this does not seem to affect the results obtained, as the impact of 
the long wavelength (lower frequency) part of the spectrum to heating is weak and as 
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the short wavelength slow waves are efficiently dissipated in the collisionless regime too, 
depositing the energy near the loop footpoints. 

Therefore, the theoretical results on the slow wave dissipation are consistent with the 
phcnomcnological results on EUV coronal loop heating and with the observed amplitude 
of the waves in the long wavelength part of the spectrum. The wide spectrum slow 
magnetoacoustic waves can, indeed, deposit sufficient amount of energy near the loop 
footpints, which is in agreement with the observations of Aschwanden et al. (1999, 2000, 
2001). We would like to emphasize that this possibility is purely theoretical due to not 
yet fully observed spectrum of the slow magnetoacoustic waves. However, our results 
indicate that the slow waves cannot be excluded from the consideration as a possible 
mechanism for heating of the coronal loops. 

Acknowledgements. The authors would like to thank Markus Aschwanden and Erwin Verwichte 
for valuable discussions. 
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